Relativistic effects and angular dependence in the reaction pp n tt"' 
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We present a new fit to the LEAR data on pp — > 7r~7r^ differential cross sections and analyzing 
powers motivated by relativistic considerations. Within a quark model describing this annihilation 
we argue, since the pions are highly energetic, that relativistic effects cannot be neglected. The 
intrinsic pion wave functions are Lorentz transformed to the center of mass frame. This change in 
quark geometry gives rise to additional angular dependence in the transition operators and results 
in a relative enhancement of higher J > 2 partial wave amplitudes. The fit to the data is improved 
significantly. 
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I. INTRODUCTION 

We recently studied effects of particle size and final- 
state interactions in the reaction pp — > tt^tt"*" within the 
framework of a constituent quark model . The aim was 
to improve previous attempts to describe the LEAR 
data on pp — > 7r"7r+ differential cross sections da/dQ and 
analyzing powers Aon in the momentum range from 360 
to 1550 MeV/c jj]. As of yet, theoretical approaches, 
whether using a baryonic or a quark picture H, |E IE 
0, 0, 0], have not been successful in reproducing the 
characteristic double-dip structure of the observables 
nor the da/dQ forward peaks at low momenta. The large 
variations of the LEAR observables as a function of the 
angle portends the presence of a substantial number of 
partial wave amplitudes already at low energies. In fact, 
the experimental data on differential cross sections as well 
as those on asymmetries point to a significant J = 2, J ~ 
3 and even higher J contributions (llj, IjJ, [ij, llj] . Model 
calculations |EIEIESS0|j however, lead invariably to 
scattering amplitudes which are strongly dominated by 
total angular momentum J = and J = 1. This is due 
to a rather short range of the annihilation mechanism in 
these models. 

In order to study possible higher partial wave J > 2 
contributions, the role of final-state tttt interactions has 
been investigated in Ref. P, H, Q, For example, in 
Ref. T] we made use of the non-planar R2 quark-flow di- 
agrams, in which a qq pair in either a '^Si or a ^Pq state is 
annihilated and momentum is transferred to a remaining 
quark or antiquark as discussed in Ref. 0, Q . Switch- 
ing on the TTTT interaction moderately improved the fit of 
in particular the double-dip structure at low mo- 
menta is slightly more pronounced. This is caused by a 
readjustment of the strengths of the helicity amplitudes 
of different total angular momentum J. Nonetheless, the 
predictions for da/dQ showed only a modest improve- 
ment over the model without final-state interactions. 
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The main improvement obtained in Ref. 1], however, 
is due to a different aspect, namely a readjustment of 
the size parameters of the intrinsic proton and pion wave 
functions, so that the radii of the proton, antiproton, 
and pions increase. In a final fit these radii are larger 
by about 7% than the respective measured charge dis- 
tribution radii. This runs contrary to the view that the 
valence quarks occupy a smaller volume than indicated 
by the charge radii. Nevertheless, introducing larger radii 
as well as final-state interactions, improves the quality of 
the fit in Ref. jlj dramatically - the forward and back- 
ward peaks of da/dfl are very well reproduced and so 
are the characteristic double-dip structures of Ao„. It 
should come as no surprise that the relative contribution 
of the higher J > 2 amplitudes in Ref. Q turns out to 
be significantly larger than in Ref. Q 

In this paper, we report an alternative approach that 
also leads to enhancement of the higher partial waves. 
It addresses the relativistic effects due to Lorentz trans- 
formed intrinsic pion wave functions in the reaction 
pp — > 7r^7r+. The reason for this rather different ap- 
proach stems from the observation that at the center- 
of-mass (cm.) energies ^/s considered in the LEAR ex- 
periment pp — !■ 7r~7r+, the produced pions are highly 
energetic. The relativistic factor in this energy range is 
7 = Ecra/'^fTi-n — 6.8 — 8.0, which means the pions are 
ejected at speeds v ~ 0.98c. One therefore expects con- 
siderable relativistic effects due to their modified internal 
structure. In the cm. frame in which the transition am- 
plitudes are calculated, the intrinsic pion wave functions 
employed must be Lorentz transformed from the pion rest 
frame. In our model the intrinsic pion as well as 

the proton and antiproton wave functions are of Gaus- 
sian form in their respective rest frames. We neglect the 
distortions of the proton and antiproton intrinsic wave 
function due to their much smaller kinetic energy and 
larger mass. 

The relativistic effects have several consequences. 
First, the intrinsic geometry of the pions is modified — 
instead of spherical particles one deals now with highly 
flattened ellipsoids in the cm. frame. Obviously, the 
reaction geometry is altered also since the boosted pion 
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wave functions enter the computation of the transition 
operators, which are obtained from an overlap integral of 
initial antiproton-proton, final pion-pion wave functions 
and the '^Si or '^Pq annihilation mechanism. Details of 
this calculation were presented in Ref. . Secondly, due 
to less overlap of the pion and antiproton-proton intrinsic 
wave functions, the annihilation range actually decreases 
but this effect is far from isotropic. As a result the tran- 
sition operators exhibit significant additional angle de- 
pendence. Finally, the annihilation amplitudes, already 
non-local in the case of non-relativistic wave functions, 
become explicitly dependent on the cm. energy ^/s via 
the boost factor 7. 

In Ref. it was shown that the relativistic trans- 
formation of the spatial part of the intrinsic pion wave 
functions introduces new angle dependent terms in the 
transition operators. These new terms also depend on 
the boost factor 7. In this paper, we present a new fit 
to the LEAR data on pp 7r~7r"'" using the R2 transi- 
tion operators of Ref. [l5j , which supports our claim that 
relativistic considerations lead to a strongly modified an- 
gular momentum content of the helicity amplitudes for 
< J < 4. Especially the amplitudes for J > 2 are 
amplified considerably. Even though we keep the origi- 
nal particle radii as in |^|^, i.e. smaller than the cor- 
responding charge distribution radii, we achieve a very 
good reproduction of the LEAR observables da/dfl and 
Aon, comparable to the one in Ref. 0. 



II. RELATIVISTIC MODIFICATIONS 

We summarize the effects of Lorentz transforming in- 
trinsic pion wave functions on pp tt^tt^ annihilation 



operators within the constituent quark model. The ac- 
tual details may be found in Ref. In the pion rest 
frame, the radial part of its intrinsic wave function is 
described by a Gaussian of the form 



■ip^(ri,rj) = A^^exp 



/3 



R.) 



(1) 



which reads in the cm. frame 



-expJ-^ E [(ra-R.)i+7'(ra-R.)f] 1(2) 



The vectors and Vj are the quark and antiquark coor- 
dinates, = 5(1"^ -I- Vj) is the pion coordinate, and (3 
is the size parameter. As mentioned above, in previous 
work U we varied the value of (3 in order to obtain an 
increase in the annihilation range. In this paper we use 
throughout the fixed value /3 = 3.23 fm~^, which corre- 
sponds to a pion radius of 0.48 fm ■ The new 
normalization factor N^^ = ^/7 A^^ is due to the condition 
that i>Tr be normalized to unity in the cm. 

The cm. wave functions of Eq. (|2Jl are used to compute 
the transition operators for the '^Pn and Si annihilation 
amplitudes. It is shown in Ref. [13 that these operators 
acquire new terms if compared with the non-relativistic 
expressions and which manifestly introduce additional 
angular dependence. We here briefly recapitulate the re- 
sults for the R2 diagrams from The complete form 
of the R2 transition operators for the vacuum ^Pq ampli- 
tude is 



f (3Po) = iA/" Aycr-R'sinh(CR-R') + Bycr • R cosh(C R-R') -f Cv(cr-R') i?cos6icosh(CR-R') 

X exp{ylR'^ + BH^ + DB? cos^ 6}. (3) 
The total "^^i amplitude for the longitudinal component is given by 

f(^Si)^iN \Al a- ■ R' sinh(C R • R') + BlO" • R cosh(C R • R' + Cl (cr • R') i? cos 61 cosh(C R- R') 



and 



exp{AR'^ + BR2 + DR2 cos^ 0} (4) 



T{^Si) =7V ATcr-R'cosh(CR-R') + BTcr.Rsinh(CR-R') +CT(cr-R')i?cos6isinh(CR-R') 

X cxp{ylR'^ + BB? + DR^ cos^ 6} (5) 



for the transversal component. The factor A/" is an overall spectively. The angle 6 is between R' and R in the cm. 
normalization and R' = R^r-— R7r+ ^md R = Rp — Rp are frame. In the experiment it is the cm. angle between the 
the relative 7r~7r"'" and antiproton-proton coordinates, re- antiproton beam direction and the outgoing 7r~ direction. 
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FIG. 1: Differential cross section and analyzing power of the reaction pp ■n~Ti'^ at T\^\^= 66.7 MeV (piab= 360 MeV/c). 
Experimental data are from Hasan et al. 4]. Solid curves denote quark model predictions with boosted intrinsic pion wave 
functions; dashed lines represent predictions with original (spherical) pion wave functions. 



The new terms mentioned previously are Cy, Cl, Ct, 
and D. The selection rules for the R2 diagrams discussed 
in Ref. Q are unchanged. More precisely, T{^Pq) and 
f{^S^) act in pp states with J" = 0+,2+,4+,... waves 
while f{^Sl) contributes only to = l-,3",5~,... 
waves. The explicit expressions for the coefficients of 
Eqs. lO"© E^re (with a being the proton size parameter 
defined similarly to (3 in Eq. (^): 



A = 

B = 

C = 

D = 

Av = 
Bv = 
Cv = 

Al = 
Bl = 



a{ha + 4(3-/'^) 
" 2(4a + 3/372) 
3(7a2 + I8a(3j^ 



8(4a + 3/372) 
3a(a + /372) 



D 



2 (4a + 3/372) 
9/3(72 - 1) 



1 



(5a + 4/3)(5a + 4/372 



a(a + /372) 
4a + 3/372 

3(a + /372)(5a + 3/372) 

2(4a + 3/372) 
3/3(72 - 1) 



-Cv 



-Av 

9(a + /372)2 
2 (4a + 3/372) 
3/3(7' - 1 



Cl 



1 - 



(5a + 4/3)(5a + 4/372) 



(5a + 4/3)(5a + 4/372) 



(6) 
(7) 
(8) 
(9) 

(10) 

(11) 
(12) 

(13) 
(14) 

(15) 



At 
Bx 



-2Av 

3a (a + /372) 



Ct — — 



4a + 3/372 
3/3a2(72 



Ct 
-1) 



(5a + 4/3)(5a + 4/372 



(16) 
(17) 

(18) 



Clearly, Cv , Cl , Ct , and D vanish in the non-relativistic 
limit 7 ^ 1 . A noticeable feature is that some of the new 
coefficients become relatively large for large 7 values. In 
particular, the DR^ cos^ 9 term dominates the exponen- 
tial part of Eqs. and, for example, Cv ^ Av,Bv 
and Cl ^ Al,Bl- On the other hand, one finds that 
Ct — At — Bt- 

At this stage we will neglect additional relativistic cor- 
rections due to the (anti)quark spinors in the final state 
pions since it was argued in Ref. [13 that boosting of spin 
wave functions results in relatively smaller corrections. 



III. RESULTS 

The LEAR data for the da/dO, and have been 
fitted using the transition operators in Eqs. ||2Jl, 101, and 
Q. As before, taking only the R2 topology into account, 
we use a distorted wave approximation (DWA) 



T' = / dRdR' 0^,(R') fR2{R, R': 



pp 



■ j±i 



(R) (19) 



to obtain the scattering matrix elements and intro- 
duce final-state interactions as in Ref. 0. The pion 
wave function (/);J^(R') is obtained from the coupled- 
channel model of Ref. 0| and can be parameterized 
with phases Sj and inelasticities rjj. The initial pp wave 
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FIG. 3: As in Fig.[T]but for Ti^b^ 219.9 MeV (piab= 679 MeV/c). 



function 4"^^ ''"'^^(R) is taken from an optical potential 

model 0. Both 0;J^(R') and *pp'='^=^^(R) still contain 
the angle dependence appropriate for the values of J and 
I. 

We select from the specific energies where LEAR data 
are available Q the same set of five energies as in Ref . Q : 
Ti^b = 66.7, 123.5, 219.9, 499.2, and 803.1 MeV corre- 
sponding to antiproton momenta of respectively piab = 
360, 497, 679, 1089, and 1467 MeV/c. Tiab is the labora- 
tory kinetic energy of the antiproton beam. In fitting the 
data, we proceed as in Ref. One starts with 7r~7r+ 
final-state plane waves, i.e. with Sj = and rjj — 1, and 
varies these parameters for < J < 4. The other param- 
eters are the relative (complex) strength A as defined by 



the R2 transition amplitude 

T;„,(R' R) = ^0 [T(3po)(R',R)-f AT (35i)(R',R)], (20) 

and the overall normalization Nq. The size parameters a 
and P of the quark model are kept fixed at the original 
values a = 2.8 fm-^ and /3 = 3.23 fni'^ ^^Q^. The 
aim is to find out how much improvement is obtained 
solely from introducing final-state interactions. 

At this initial stage, no relativistic effects have been 
included which implies 7=1. The results of this fit are 
plotted in Figs. ^0 as dashed lines. They will serve as 
reference point and were obtained previously in Ref. 
where the model with and without final-state interactions 
was discussed. It was found that the improvement from 
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FIG. 5: As in Fig.Elbut for Ti^^= 803.1 MeV (piab= 1467 MeV/c). 



just final-state interactions is very modest — the double- 
dip structure in the analyzing power is somewhat more 
enhanced but the forward peaks in the differential cross 
sections at Tiab = 66.7 and 123.5 MeV (piab = 360 and 
497 MeV/c) as well as the backward peaks at Tiab = 499.2 
and 803.1 MeV (piab = 1089 and 1467 MeV/c) are poorly 
reproduced. 

Next, we consider the effects of relativistic distortions 
of the intrinsic pion wave functions. Hence, the tran- 
sition operators of Eqs. ®-© must be used. The ex- 
ponential part is dominated by the D term (typically 
D ~ 200 fm~^, about two orders of magnitude larger 
than the C term), which makes the partial wave am- 
plitudes T"^ very sensitive to the value of the relative 



strength A = \X\exp{i6\) of the versus the '^Pq 
mechanism and to the pion wave parameters 77 j and Sj. 
This is a consequence of the strong angle dependence of 
exp{DR^ cos^ 6}. In order for the partial wave expansion 
of the scattering amplitude T to converge, all ingredients 
must be fine-tuned, in particular the phases Sj and in- 
elasticities rjj for < J < 4. It should be emphasized 
that a and (3 remain fixed. 

The improved new fit is illustrated in Figs.ni^l where 
the predictions that include relativistic effects and final- 
state interactions are presented as solid curves. The new 
fit is a significant improvement over the non-relativistic 
version of the R2 annihilation model. The differential 
cross sections da/dQ are particularly well reproduced 
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Tlab [MeV] 


66.7 


123.5 


219.9 


499.2 


803.1 


Plab [MeV/ 




369 


497 


679 


1089 


1467 


7 = y/s/2m 




6.841 


6.940 


7.106 


7.564 


8.033 


|A| 


1.165 


0.889 


1.090 


1.377 


1.320 




154.51 


149.13 


82.57 


184.14 


184.76 


iVo [10' 1 


2.999 


3.061 


3.068 


1.739 


4.999 


TABLE I: Quark model parameters as 


function of Tj^b- 


Tlab [MeV] 


66.7 


123.5 


219.9 


499.2 


803.1 


[MeV] 


1910. 


1937. 


1983. 


2111. 


2242. 


no 


1.00 


0.716 


0.811 


1.00 


1.00 


So 


67.00 


54.37 


41.42 


-1.704 


-17.15 


Vi 


1.00 


0.583 


0.666 


0.999 


1.00 


Si 


43.68 


37.11 


42.62 


-31.62 


21.56 


m 


0.968 


0.973 


0.992 


0.538 


0.624 


&2 


-17.06 


-16.55 


-18.78 


-37.34 


-13.51 


V3 


1.00 


0.923 


1.00 


0.946 


1.00 


Sa 


-11.93 


-13.80 


-16.07 


-6.873 


-11.91 


Vi 


1.00 


1.00 


1.00 


0.946 


1.00 


Si 


-0.03 


-0.03 


-0.04 


-0.038 


-0.07 



TABLE II: Phases shifts Sj and inelasticities i]j of the final- 
state TTTT interaction for < J < 4 as function of ^/s. 

for the energies Tiab = 66.7, 123.5, and 219.9 MeV 
(piab = 360, 497, and 679 MeV/c). For Tiab = 499.2 and 

803.1 MeV (piab = 1089 and 1467 MeV/c) the results 
for da/dn show improvement in forward and backward 
direction. At the two higher energies, the experimental 
double-well structure of the cross section is reproduced 
but some problems remain near 9 = 7r/2. 

As for the analyzing powers Aq^, we are able to repro- 
duce their double-dip structure, already present in the 
LEAR data at the lowest energies. When the experimen- 
tal errors are small, like for instance for Tiab = 219.9 and 

499.2 MeV, the fit is particularly successful. At lowest 
energy, Tjab = 66.7 MeV, the error bars in Aon are rather 
large and we do not fit very well certain data points. 
Overall the improvement is clearly important and of sim- 



ilar quality as the results of Ref. [l| which were obtained 
from increasing the particle radii. 

In Tab. HI we give the values of the quark model pa- 
rameters as a function of Tiab while the tttt phases Sj 
and inelasticities r]j are conveniently listed as a function 
of cm. energy ^/s in Tab. ^ In the quark model pa- 
rameters of Tab.Uone notes that the relative strength |A| 
of the ■^5*1 versus the ^Pq mechanism varies with energy 
and the average value |A| ~ 1.17. There is no a priori 
reason why A should be energy independent. The energy 
variation of the inelasticities rjj in the tttt interaction ob- 
served in Tab.HTlmav be associated with the presence of 
TTTT resonances, in particular with the /2 at 2010 MeV in 
the J — 2 amplitude. 



IV. SUMMARY AND CONCLUSION 

We present a new fit to the LEAR pp 7r^7r+ data 
based on relativistic effects in the annihilation. These ef- 
fects are due to Lorentz transformations of intrinsic pion 
wave functions from the pion rest frame to the cm. frame 
of the reaction. The impact from the modified geometry 
of the pions is considerable and carries with it a much 
richer angle dependence of the transition operators. 

The quality of this new fit is very similar to the one 
of a previous fit Jj. However, the approaches are very 
different. While in Ref. we achieved a good fit by 
means of particle radii increase, so the antiquarks and 
quarks do overlap more, in this work we keep the original 
particle sizes fixed but take into account the relativistic 
distortions of the pions in the cm. frame of the reaction. 

We conclude that even though observed data at LEAR 
on da/dfl and Aon indicate strongly the presence of 
higher partial wave amplitudes and therefore seem to call 
for an increase of the annihilation range, one can obtain 
a similar improvement by additional angle dependence 
in the transition operators. This additional angle depen- 
dence arises naturally when boosted intrinsic pion wave 
functions are employed in the quark model calculations, 
as becomes clear from the present fit. Relativistic ef- 
fects in the reaction pp 7r^7r+ therefore should not be 
ignored. 
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